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We consider homogeneous non-abelian vector fields with general potential terms in an expanding 
universe. We find a mechanical analogy with a system of N interacting particles (with N the 
dimension of the gauge group) moving in three dimensions under the action of a central potential. 
In the case of bounded and rapid evolution compared to the rate of expansion, we show by making 
use of a generalization of the virial theorem that for arbitrary potential and polarization pattern, the 
average energy-momentum tensor is always diagonal and isotropic despite the intrinsic anisotropic 
evolution of the vector field. We consider also the case in which a gauge-fixing term is introduced in 
the action and show that the average equation of state does not depend on such a term. Finally, we 
extend the results to arbitrary background geometries and show that the average energy-momentum 
tensor of a rapidly evolving Yang-Mills fields is always isotropic and has the perfect fluid form for 
any locally inertial observer. 



I. INTRODUCTION 

The possibility that some of the unknown components 
of the universe at different epochs (inflaton field, dark 
matter or dark energy) could be described by means 
of homogeneous vector fields rather than scalar fields 
have conflicted traditionally with the stringent limits on 
isotropy imposed by CMB observations [THH]. The only 
scenarios in which these limitations could be evaded are 
those in which either only the temporal components of 
the vector fields are present or, in the case in which spa- 
tial components are also evolving, if the vector configu- 
ration guarantees an isotropic energy-momentum tensor. 
Two possible isotropic configurations discussed in the lit- 
erature are the presence of a triad of mutually orthogonal 
vectors [TUHH] or the existence of a large N number of 
randomly oriented vectors, so that the average isotropy 
violation is kept relatively small of order 1 / y/~N [5] . 

However, recently [13] a new possibility was devised 
in which a general isotropy theorem for vector fields was 
proved. Provided vector field evolution is bounded and 
rapid compared to the rate of expansion, it was shown 
that the average energy-momentum tensor is always diag- 
onal and isotropic for any kind of initial configuration. A 
typical example of rapid bounded evolution corresponds 
to massive vector fields with masses larger than the Hub- 
ble parameter. The theorem was proved in the case of 
abelian fields with standard Maxwell kinetic term. In this 
work we are interested in extending it to more general 
non-abelian gauge theories in which gauge-fixing terms 
could also be present. 

Unlike abelian theories, the presence of multiple vector 
fields in Yang-Mills theories allows to implement in a nat- 
ural way the triad configuration mentioned before. Thus, 
in [TTJ [T2] it was found that for SU(2) groups, homoge- 
neous and isotropic configurations compatible with FRW 
backgrounds can be parametrized with a single function 
of time. Different proposals for dark energy or inflation- 
ary models based on this triad solution with different 
types of actions have been considered in the literature 
[5] QTJ [T2J Q3] . However apart from this particular type 



of solution no further isotropic configuration has been 
proposed to date. The theorem presented in this work 
applies quite generally without resorting to particular so- 
lutions by means of a generalization of the classical virial 
theorem and allows to ensure homogeneity and isotropy 
in average irrespective of the concrete dynamics of the 
non-abelian fields for arbitrary semi-simple groups. 

The paper is organized as follows, in Section II we in- 
troduce Yang-Mills theories in FRW backgrounds and ob- 
tain the corresponding equations of motion and energy- 
momentum tensor components. In Section III we intro- 
duce the generalized version of the virial theorem and 
show that in average the energy-momentum tensor is di- 
agonal and isotropic for semi-simple groups. As an exam- 
ple we compute the average equation of state for certain 
SU(2) solutions. In Section IV we extend these results 
to actions including gauge-fixing terms and obtain the 
average equation of state in different cases. We also con- 
sider the case with temporal components only. Section V 
is devoted to the extension to more general background 
metrics and finally Section VI contains the main conclu- 
sions of the work. 



II. YANG-MILLS THEORIES 

Starting from a theory whose action has a global sym- 
metry, if we would like to make it invariant under local 
symmetry transformations, we must add new bosons and 
define a covariant derivative: 

= d M +A^; A^ = -igA a J a , (1) 

where T a £ Q , a = 1 . . . N are the symmetry group 
generators. We will assume a compact semi-simple Lie 
group with a finite number of generators. In this case, 
it is always possible to find an orthonormal basis for the 
group generators for which: 

Tr(T a T b ) = ^S ab . (2) 
and the structure constants defined as: 
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[T a ,T h ] =ic abc T c 



(3) 



are totally antisymmetric (see |15|). This property will 
be very helpful for our future discussion. 

In order to write a kinetic term for gauge bosons, we 
consider the curvature tensor associated to the covariant 
derivative: 



F MV = [£>„, D v ) = d^A v - duA^ + [Aft, A,} 



(4) 



Notice that this expression is still valid in curved space- 
time due to the antisymmetry properties of the ten- 
sor. The corresponding components read: 



expected for the standard kinetic term we are using. On 
the other hand, the energy-momentum tensor is given by: 

T% = (^F^F a af3 + V^j S f " l/ -F a ^ a F^ a ~2V / M ab A a ^A b l/ . 

(12) 

Expressed in components, 



j— — ■ {AfA'l; + 2gc a bcA1 A\A\ 
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(13) 



F„ 



= d^-d^+gCabcAlAl. 
Finally, the Yang-Mills lagrangian density reads: 



(5) 



-^kinetic r\ o ^-^uv-^^ ) 

2g l 
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In this work we will study Yang- Mills theories with a 
potential of the form V(M a bA a p A bp ) 1 with M a b a constant 
symmetric matrix, so that gauge symmetry is explicitly 
broken. Thus, the action in general curved space-times 
will read: 

S = j ' d'x^g (^-- A F« u F a ^ -V{M ab A a p A b ^ . (7) 
The corresponding equations of motion are given by: 



Pk 



Tk — „ a i \ {MM + ^gcabcAf A b Q Al 



2a 4 (77) 
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^fc^fc + ^gCabcA^A^AQ 



a 4 (?l) 

+ g 1 CabcCade (A\A* - A)A$ A%A%) - V 

A a A a 

- 2V'Mab^r^ , k= 1, 2, 3; 



a 4 ( V ) 



AfAj — gcabc (AfA b + AjA b ) Aq 



(14) 



+ g 2 CabcCade A b A jAqAq — g 2 CabcCadeM^kM A% 



+ 2V'Ma b a 2 {r,)A a l A b A , i ? j 



j k 
(15) 



pa \v 



gCabcFl v A c v + 2V'M ab Al = , 



(8) 



dV{oi 



dx 



where V 

We will be interested in cosmological solutions for 
which the gauge fields will depend only on time (77) 
and the metric tensor will be given by the that of a flat 
FLRW metric: 



ds 2 = a 2 (77) (drj 2 -dx 2 ) 



(9) 



rpQ _ _ 



A'f 



+ 2V / M ab a 2 (v)A b 



gCabc-AjA^ + g 2 CabcCbdeAftA'j Aj 



(16) 



The average of the pressure along the three spatial direc- 
tions will be denoted by: 



(17) 



where 77 denotes the conformal time coordinate. The field 
equations expressed in components can be written as 

gCabcA\A\ + g 2 CabcCbdeA*A\A\ + 2V'M ab a 2 { V )A b =0, 

(10) 

for n = 0, and 

M - 9Cabc (2A b A c + A b A^j + g 2 cabec b de (A$ : A e A c 
- A d i A e j Af ) -2V'M ab a 2 {ri)A\ = Q, (11) 

for (i = i. Notice that there is no second time deriva- 
tive of the temporal component in these equations, as 



It is possible to show that, using the equations of mo- 
tion of the temporal components (10), the energy fluxes 
vanish, indeed: 



T? = 



(18) 



III. GENERALIZED VIRIAL THEOREM FOR 
NON-ABELIAN FIELDS 

We will follow the method presented in [13], in or- 
der to compute the average components of the energy- 
momentum tensor. With that purpose, we will general- 
ize the virial theorem in order to apply it to non-abelian 
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fields, thus we define: 

Aa Ab 
r ab = A i A j 

ij a 4 (77) 



i,j = 1,2,3; a,b= 1...N . (19) 



Now Af is alike the "i" position coordinate of the 
classical point particle "a". Thus, in the case Aq = 0, 
the problem is analogue to a mechanical system of N 
interacting particles in three dimensions, with N the 
group dimension. 

Assuming a rapid evolution of the vector field, we 
can neglect the expansion of the universe at the field 
equations, the time derivative of the previous expression 
becomes: 



Aa Ab 
Qab = A i A j 



A? A) 



lJ a 4 ( V ) a 4 (r,)' 



hi = 1,2,3; 



a, b = 1...N. (20) 



Integrating the last expression on [0, T] where T is larger 
than the typical time scale of the vector field evolution 
lo^ 1 but smaller than the typical time scale of the uni- 
verse expansion H~ , i.e. H^ 1 »T» w -1 , we obtain: 



Gf b (T) - G°"(0) = I A- A) A* A] 
T \ a 4 (77) a 4 (77) 



(21) 



If the motion is periodic or it is bounded, the left hand 
side of the equation vanishes and therefore: 



Aa Ab AO. Ab 

i 3 , * 3 



. 



(22) 



Using these equations, we will show in the following that 
the average energy-momentum tensor is diagonal and 
isotropic. 



From the off-diagonal part of the tensor (151, noting 
that: 



A(i A< j) ~ '^9 c abcA^ i A^AQ — g 2 C abc Cbde (^A^A^Aq 



A' n 



- A%A^AIA%) + 2V'M ab a 2 (r 1 )A\ i A a j) , (23) 

where the parenthesis in the sub-index means sym- 
metrization, we can write: 



^ = Jw) + ***** 



A a A a -4- A a A a 
a 4 (ry) \ A « A f) +A V A i) 



with i =/= j 



(24) 



Using the average expression (22), we directly obtain: 
1 



A a A a 



A a A a 
A (i A 3) 



0, i+j. (25) 



We can prove the isotropy of the diagonal elements in a 
similar way. Indeed, the diagonal spatial components of 
the energy-momentum tensor are 



1 

4^ 



pa pa pX + y 



1,2,3; (26) 



where t\ is the anisotropic part of these components, 
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A kMi — ^9 c abcAkA b k AQ 
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(27) 



which is equal to the expression Tj when i = j. Therefore 
1 



' 1 



a 4 ^) 



A a A a 
A (k A k) 



A tk A l) 



(4)=0, k = 1,2, 308) 



where again we have made use of (22). Thus, we can 
conclude that the virial theorem guarantees the isotropy 
of the energy-momentum tensor and its diagonal form in 
average. 

A. Equation of state for SU(2) particular solutions 

For SU(2) groups some particular solutions are already 
known in flat space-time, see [10] where a triad solution 
was proposed: 



At = Km 

A a Q = . 



with k a -k b = S ab 



(29) 
(30) 



Notice that in those solutions the temporal components 
of the vector fields vanish. For abelian field theories the 
vanishing of the temporal component comes naturally 
from the equation of motion, as it can be seen in ( 10 ) 
by making g = 0. In the non-abelian case it would be 
possible to have solutions with non-vanishing time com- 
ponents, however in the present example we are not con- 
sidering this possibility. Notice also that this kind of 
solutions only depend on a single function of time. We 
have checked that for the above ansatz, solutions also 
exists in an expanding background. 
Let us define: 

V = ^ )Cabc c ade (A b Ai){A^) + V (-M ab ^§ 

( 31 ) 

The average energy-momentum tensor components then 
take the form: 



(p) 




(Pk) 


= (p) 


pX) 


= 0, 


(if) 


= . 



1 A? A? 



lAp^l 
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V 



(32) 

(33) 

(34) 
(35) 
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For simplicity we will assume a power law potential: 

V (M ab A a p A b p ) = - {-M 2 A a p A ap ) n . (36) 



For the particular ansatz (29), we get: 

'3 / 2 , (3M 2 / 2 )" , 3 g 2 



(p) 



(P) 



2a 4 (ri) 2a 2n (r]) 2« 4 (i)) 



3 f 



(3M 2 / 



2 r2\n 



3 g 2 



2o 4 (»)) 2a 2n (?7) 2a 4 (ry) 



27 = 0, 



/ 4 



/ 4 



We can write: 



3/ 2 



(Tj) = 



{P+P) = (l + lp)(p), 



(37) 



(38) 



(39) 
(40) 



(41) 



where 7 is the average value over an oscillation and j p 
the periodic part of the oscillation, that can be neglected 
(see [H]). Calling r = *- and taking the scale factor as a 

constant r = we can define an effective potential: 



y e// (r) = i(3MVr + ^V. 



such that: 



P 



2 a 2 (77) 



Veff(r) • 



(43) 



Notice that the non-abelian contributions appear as a 
quartic term in the effective potential, so that in the 
strong coupling limit, we expect the equation of state 
to correspond to that of a scalar field with a quartic po- 
tential, i.e. that of radiation [T5] . 

In order to compute the equation of state, we first ob- 
tain 7 following [Tni as: 



2 Sll dry/p-V eff (r) 

Jr ~ y/p-V Bff (r) 



(44) 



where r + and r_ are the turning points of the effective 
potential where r = 0. The corresponding equation of 
state is (p) = w(p), being w — 7 — 1. The integrals 
cannot be computed analytically but the numerical 
results 1 show that w depends on n, p, g and M varying 



1 Now, unlike the abelian case, fdV e ff/df is not proportional to 
V e ff, except when n = 2, and that is the reason why the virial 
theorem can not be used to obtain the equation of state analyt- 
ically as done in 1131 




Log 10 (£) 

FIG. 1: In this figure, it is shown how the equation of state 
(p) — u)(p) changes by varying g in the SU(2) non-abelian case 
with potential ( 36 1 , for a given energy density and M param- 
eter. The different lines correspond to different values of the 
potential exponent n in ( 42 1 . We can observe that for large 



g values, the non-abelian term dominates and the equation 
of state approaches the radiation behaviour independently of 
n. On the contrary, for small g, the abelian case is recovered 



(42) between radiation for large g (the new non-abelian 



term dominates in (42)) and the corresponding uj of the 



abelian case |13) for a given n when g is small (see Fig. 
0. 

The behaviour of the equation of state depends 
also on the energy density. If the energy density 
decreases as the universe expands so that the amplitude 
of the field oscillations decreases in time, then at late 



times the smallest power in the effective potential (42) 
will dominate. Thus, for power-law potentials we can 
have two possibilities: for n > 2, the potential term will 
contribute at high energies (early times) and ui = ^=j, 
whereas at low energies (late times), the non-abelian 
quartic term will dominate and uj = |. For n < 2, 
the situation is reversed and at high energies, we have 
domination of the non-abelian term and w = | , whereas 
at low energies uj = 2=j (see Fig. [2] and reference [TB"]). 



IV. YANG-MILLS THEORY WITH A 
GAUGE-FIXING TERM 

As is well-known the standard kinetic term for gauge 
fields ^ does not contain a time derivative of the tempo- 
ral components which implies that no conjugate momenta 
can be defined for those fields. Accordingly, temporal 
components cannot be quantized along with the spatial 
ones, thus preventing an explicitly covariant quantiza- 
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FIG. 2: In this figure, it is shown how the equation of state 
(p) = u)(p) changes by varying p [M 4 ] in the SU(2) non- 
abelian case, for g = 50, n = 1 and M = 1. It can be seen how 
the behaviour of the higher power dominates at high energy 
density while the contrary happens when the energy density 
is low. 



tion of the theory. Covariant quantization thus requires 
the modification of the kinetic term by including new 
gauge-breaking terms |17j . Those terms can be included 
directly in the action as in the Gupta-Bleuler formalism 
or they could appear as a consequence of the selection 
of a representative within each gauge orbit in the path 
integral approach. Restricting to the case of quadratic 
gauge-breaking terms the action reads: 

C = -\f%F° + ^-{V p A a P f - V(M ab A a p A b ») . (45) 

We can express the covariant equations of motion as 

F^-gc abc F b u A c "+2V'M ab A b l H {A a p % = • ( 46 ) 

By introducing the homogeneity condition, this expres- 
sion is reduced to: 



fig + gcabcMAi + g 2 c abc c bde A d A e i A c i 



(47) 



+ 



a a z 



5 ab + 2V , M ab a 2 (r,))A b =0 1 



A? 



gc abc [2AlAl+AlAl 



A a A e 4 CX 
A i A 3 A 3J 



' ' ' 9 2c abc c bde (Af AqAq 

2V'M ab a 2 (r,)A\ = . (48) 



In order to have oscillations we must restrict the param- 
eter space and initial conditions. In the case of a power 
law potential like (36) with an odd power, those condi- 



tions simply require f < 0. Notice that we do not need 
to impose any condition over the temporal component in 
order to ensure isotropy. 



The energy-momentum tensor has the form: 
1 



pa pa p\ + y 



pa ppa 



p\ r f v J 9p,f * fj, * up 

2V'M ab A%Al + | [S^ [{V p A a P f 
2A a x W x (W p A a P)] -4A a iti y u) (V p A a P) 



(49) 



where the parenthesis in the sub-index means sym- 
metrization. If ^ < 0, the energy density associated to 
the temporal component of the field is generally negative 
defined. This is consistent with the standard interpreta- 
tion of Aq as a ghost field. 

In the following, we will show that the introduction 
of the gauge breaking term does not spoil the averaged 
isotropy of the energy-momentum tensor. Let us first 
consider the energy fluxes. We see that: 



Af r 
a 4 

2* 



£Aq + gC abc AjAj + g 2 Cabc c bdeAQAjAj 



S ab + 2V'M ab a 2 (r)))A b Q 



= 0, 
(50) 



where in the last step we have used the equations of mo- 
tion of the temporal components (47). On the other 



hand, the new £ term does not contribute either to T? 
or t\. Therefore, following the same procedure of the 



previous section (see Eqs. (19) to (28)), we can conclude 



that the virial theorem also guarantees in this case the 
isotropy of the energy-momentum tensor and its diagonal 
form in average if the spatial fields oscillate fast enough 
compared with the expansion rate of the universe. 



A. Equation of state 

In order to obtain the equation of state we will proceed 
in a different way as it has been done in the previous sec- 
tion. We take the trace of the energy-momentum tensor: 



rpf, 
P 



p-3p = W- 2V'M ab A a ^ 

2£ [(V p A ap ) 2 + A a ^V ^{V p A a p )} . (51) 



with p the average pressure along the three spatial direc- 
tions defined in (17). Thus, we have: 



3p = W - 2V'M ab A a 11 A\ 
+ 2£V A (A a x W p A a p) . 



(52) 



Notice that the fast oscillations condition is not manda- 
tory for the temporal part since it does not contribute 
to the anisotropics. Nevertheless, if the oscillations are 
indeed faster than the expansion rate of the universe, 
taking the average, we get: 



(i-M(p) 



(4V - 2V'M ab A a p Al) 
2£(V x (A a x V p A a P)) . (53) 



() 



Notice that the £ term, when neglecting the scale factor 
derivatives, is nothing but the temporal derivative of the 
Gqq function used in the virial theorem which in average 
vanishes, i.e.: 



y^iab 
^00 



Aa Ab 



,6 = 1. 



a 2 (77) 

and as in the spatial case (J22|), we get: 



■ N, 



(54) 



(G 



^ab\ 
00/ 



A a A b 
^O^O 



Aa \b 
^0^0 



a 4 (v) a 4 (77) 



, a, b = 1 ... N. (55) 



Thus, we can give an expression for u depending only 
on the average of the potential and the energy density: 

_ 1 (2V'M ab A a »A\ - 4V) 
W ~3 + m ' 



Considering a power law potential, the expression reduces 
to: 



2(n - 2) (V) 
3 (P) 



(57) 



As it can be seen the solution for a potential with n = 2 
can be obtained analytically, it behaves as radiation uj = 
I . That is the expected result considering the behaviour 
of the same case without gauge-fixing term, see Fig. [T] 



Introducing this expression in (57) 
result 



as [161 

theory without gauge-fixing term 



we reach the same 
for a scalar field and [13] for an abelian 



1 



(63) 



2. Non-abelian case 

For non-abelian theories we must compute the aver- 
age of the potential and the energy density in (57). As 



commented before, we cannot apply the virial theorem 
as in the abelian case for simple power law potentials be- 
cause now some non-abelian terms (those multiplied by 
g 2 ) act as additional potential terms. Thus the effec- 
(56) tive potential becomes a sum of power laws, and V' e ^A 2 
is not proportional to V e ff. In any case, if we assume 
< (V) /(p) < 1, we can constrain the value of uj: 



2n-3 

3 



I , for n < 2 ; 
- ' < (64) 
^ , for n > 2 . 



B. Particular case: A£ = A%(i])5l 



1. Abelian case 

In the abelian case (or g — > 0) we can obtain the exact 
expression for the average equation of state by following 
the virial method despite Aq ^ 0. Using (55) and (22) 



and the equations of motion ( 47 ) and ( 48 ) , we have 



^0^0 



2Ma f M {V'AUl) , 



A?A? 



: -2M ab a 2 ( V )(V'A^At) , 
and the average energy density can be written as 



(58) 
(59) 



(P) 



2 a 4 (ri) 2 a 4 (r]) 
1 



V 



a\r,) 



(ZA a Q + 2V'M ab a 2 (r 1 )A b Q )A a 



(60) 



Using (47) and neglecting time derivatives of the scale 



factor, we see that the last term vanishes and finally we 
get: 



(p) = 



1A?A? ^Ag^g 
2 a 4 (rj) 2 a 4 (77) 



V 



Then using (58) and (59), we obtain: 



(p) = (n + l)(V) . 



(61) 



(62) 



Let us consider the case in which the only relevant 
component of the homogeneous Yang-Mills theory with 
a gauge-fixing term ( i.e. £ ^ ) is the temporal part of 
the fields, 



The equations of motion become: 



(65) 



with no contribution from the non-abelian terms. This 
equation is similar to the standard scalar field equation, 
although with variable mass depending on the scale fac- 
tor derivatives. However the non-zero components of the 
energy- momentum tensor are completely different, 

p = V+Uv p A?f-2V'M a bA«Al-tA a0 V {V p A a P) 



—0 , by c.o.m. 



r + | (v P A p ) 2 , 



(67) 



P = -V-^(V p AP) 2 -ZA a0 V (V p A a P) 



Adding them: 

p + p = -iA a °V (V p A a p) = 2V'M ab A a °A b , (69) 
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where in the last step the equations of motion have been 
used. We can give an expression for the equation of state: 



2V'M ab A a °A b 

uj = " — 1 

P 



(70) 



Now the oscillatory behaviour is not required as the 
isotropy is guaranteed, so for a simple power-law poten- 
tial, we get: 

2n , . 

(71) 



1 



1 . 



We see that contrary to the scalar field standard be- 
haviour, when the potential dominates u) = 2n — 1; 
whereas lj = — 1 when the potential term is negligible. 
When the field is rapidly oscillating, we can use the virial 
equations ( 63 ) and the average uj behaves as a standard 
scalar field ui 



77+1 



since 



V 



2n 



(72) 



as it can be straightforwardly deduced from Eq. ( 58 ) 



This does not mean that the temporal component be- 
comes a standard scalar field, but that in average its 
equation of state behaviour is the same. A good way to 
illustrate this is by studying the case Va = M (A^A 0, ^) 
and its equivalent V$ — Mtp 4 for a scalar field. The con- 
ditions to have oscillations are M > and £ > 0. In this 
particular case, the pressure behaves completely differ- 
ent when the field is a scalar, a temporal component or a 
spatial component of a vector field; although, in average 
they behave in the same way. 

In Figure [3J the different cases are shown for g = in 
a radiation-dominated universe. For a scalar field, the 
pressure will oscillate in the interval — p < p < p, so 
pressure is modulated by energy density. If we have a 
spatial component of a vector field (A% = 0), the average 
pressure over the three spatial directions p always equals 
p/3, and it does not oscillate. This odd behaviour can 
be explained by realising that the trace of the energy- 
momentum tensor (51) vanishes for this potential if the 
temporal component is zero, then p — 3p = must be 
satisfied at any time. Finally for a temporal component 
of a vector field (Af = 0) , pressure oscillations exceed the 
energy density value on the contrary to the scalar case. 



V. ISOTROPY THEOREM FOR YANG-MILLS 
THEORIES IN A GENERAL BACKGROUND 
METRIC 

The previous results can be directly extended to 
general space-time geometries (not necessarily homoge- 
neous) as was shown in |13j . Thus, let us consider a 
locally inertial observer at x^ = and write the metric 
tensor using Riemann normal coordinates around Xq |18j : 

9^{x) = ri^v + \R^ a upx a x 13 + . . . (73) 



Let assume that the following conditions hold: 

1. The Lagrangian is restricted to the Yang-Mills form 
with or without a gauge-fixing term. 

2. The vector field evolves rapidly: 

I^J « K) 2 , and \d,A1\ « \Af\, 

for i, j = 1, 2, 3 ; and a = 1 . . . N , (74) 

for any component of the Riemann tensor. w° is 
the characteristic frequency of Af 

3. Af and A\ remain bounded in the evolution. 

The second condition implies that if we are only inter- 
ested in time scales of order (wf) _1 , then we are in a 
normal neighbourhood and we can neglect the second 
term in ( 73 ) and also work with a homogeneous vector 



field. In such a case, it is possible to rewrite all the above 
equations with a(rj) = 1. Thus, by using an interval [0, T] 
that verifies the condition: 



I^J«T- 2 ««) 2 



(75) 



for any of the components of the Riemann tensor and of 



the vector, it is possible to obtain (21 ) and prove that the 



mean value of the energy-momentum tensor is isotropic. 
Thus, if oscillations are fast compared to the curvature 
scale, the average energy-momentum tensor takes the 
perfect fluid form for any locally inertial observer. 



VI. CONCLUSIONS 

In this work we have considered homogeneous non- 
abelian vector fields in an expanding universe with arbi- 
trary potential term. By means of a generalized version 
of the virial theorem, we have shown that for bounded 
and rapid evolution compared to the rate of expansion, 
the average energy-momentum tensor is isotropic for any 
kind of initial configuration of the fields. This result holds 
irrespective of the complicated dynamics of the coupled 
fields in non-abelian theories and also in the presence of 
gauge-fixing terms. This opens up the possibility of using 
non-abelian fields in cosmology beyond the simple triad 
configurations considered so far. The results also hold for 
locally inertial observers in arbitrary geometries provided 
the field evolution is sufficiently rapid. 
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FIG. 3: In this figure the average pressure along the three spatial directions p (grey line), energy density (blue dashed line) 
and average pressure (p) = (p)/3 (black line) of three different cases are shown. These cases are (from left to right) a scalar field 
with potential V</, = Mcf> 4 , a vector with spatial components only and a vector with temporal component only with potential 
V A = M (A"A a ") , where M = 100, £ = 1, g = 0. The y-axis is normalized to the initial value p{t'o) and time, t', is in H Q 
units. We are considering a radiation dominated-universe. Notice that in the spatial components case (p) = p. 
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